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Abstract 

We analyze the effects of a distortion of the nuclear potential of a molecular quantum dot (QD), 
as well as a shift of its equilibrium position, on nonequilibrium- vibration-assisted tunneling through 
the QD with a single level (e^) coupled to the vibrational mode. For this purpose, we derive an 
explicit analytical expression for the Franck-Condon (FC) factor for a displaced-distorted oscillator 
surface of the molecule and establish rate equations in the joint electron-phonon representation to 
examine the current-voltage characteristics and zero-frequency shot noise, and skewness as well. 
Our numerical analyses shows that the distortion has two important effects. The first one is that 
it breaks the symmetry between the excitation spectra of the charge states, leading to asymmetric 
tunneling properties with respect to > and Ed < 0. Secondly, distortion (frequency change 
of the oscillator) significantly changes the voltage-activated cascaded transition mechanism, and 
consequently gives rise to a different nonequilibrium vibrational distribution from that of the case 
without distortion. Taken in conjunction with strongly modified FC factors due to distortion, this 
results in some new transport features: the appearance of strong NDC even for a single-level QD 
with symmetric tunnel couplings; a giant Fano factor even for a molecule with an extremely weak 
electron-phonon interaction; and enhanced skewness that can have a large negative value under 
certain conditions. 

PACS numbers: 85.65.+h, 71.38.-k, 73.23.Hk, 73.63.Kv 
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I. INTRODUCTION 



Recently, the effect of unequilibrated quantized molecular vibrational modes on 
electronic tunneling has become an active issue, both experiment allyii2i»^i^i^'^'^ and 
theoretically.-'^i^iiiii^ii^i^ii^'i^>i^ii^>i2ii22i2i In particular, electronic transport measurements 
in suspended Carbon nanotubes have yielded I-V curves demonstrating perfect signatures 
of phonon-mediated tunneling, e.g. stepwise structures with equal widths in voltage and 
gradual height reduction by the Franck-Condon (FC) factor.-i^ 

To date, most theoretical works have taken into account only the displacement of the 
molecular oscillation arising from the entrance of an additional electron into the molecule. 
It is well known in molecular spectroscopy that upon an electronic transition, the poten- 
tial surface involving the normal coordinates of the molecular oscillator generally undergoes 
a displacement, distortion, and rotation simultaneously. Accordingly, the molecular ab- 
sorption and fluorescence line shapes depend on the square of the overlap integral between 
the initial wave function of the harmonic oscillator and the displaced-distorted-rotated one 
within the Born-Oppenheimer approximation, which is referred to as the FC factor in the 
literature.— Likewise, in the event of electron tunneling, an excess electron entering into a 
molecule from electrodes can induce a shift of the equilibrium position of the internal vibra- 
tional mode (displacement of the harmonic potential surface) and a simultaneous change of 
the vibrational frequency (distortion of the potential), both of which can modify the elec- 
tronic tunneling rate and thus the nonlinear current- volt age characteristic of the molecule. 
(The rotation of the potential surface is ignored in this one-dimensional model.) J. Koch 
and F. von Oppen have studied this issue by calculating the overlap integral numerically 
in coordinate representation.— M.R. Wegewijs and K.C. Nowack have also analyzed this 
effect on electron tunneling by developing an explicit expression for the overlap integral in 
terms of Hermite polynomials.— However, this issue still requires further investigation. The 
main difficulty involves the derivation of an analytical expression for the vibrational overlap 
integral, which can be traced back to 19302^ and has generated many studies in quantum 
chemistry.— In this paper, we will employ Fan's method of integration within an or- 
dered product of operators (IWOP) to derive an explicit analytical expression for the FC 
factors between two adiabatically displaced-distorted potential surfaces in terms of gener- 
alized Laguerre polynomials.^^ The resulting exact FC factors will then be incorporated 
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with the rate equation in a joint electron-phonon representation to analyze the effects of 
distortion on sequential tunneling through a single molecule. 

Furthermore, it is reasonable to expect that this distortion effect has an even stronger 
impact on current fluctuations, which has thus far not been studied despite much re- 
cent work investigating the unequilibrated phonon effect on current noise through a single 
molecule.— ii^ii^>i^>22 Recent developments in experimental measurement techniques for meso- 
scopic electron transport make it possible to detect the probability distribution of charge 
transmitted in a fixed time interval, namely the counting statistics of fluctuating electronic 
current.— '^'^'^'^'^'^ Very recently, the full counting statistics (FCS) of tunneling current 
though a quantum dot (QD) were investigated in the Coulomb blockade regime by means of 
a rate equation.— i^^^^ Moreover, many studies have been devoted to the FCS of a single-level 
QD coupled to an oscillator with a harmonic potential.^- Therefore, another purpose of this 
paper is to study the counting statistics (up to the third moment) of electronic tunneling 
in a molecular QD based on number-resolved rate equations,— focusing on the displaced- 
distorted effect (i.e. an oscillator with an anharmonic potential), the external- voltage-driven 
unequilibrated phonon effect, and finite phonon dissipation. 

The rest of the paper is arranged as follows. In Sec. II, we describe the model Hamiltonian 
and theoretical formulation of the rate equation. We briefly discuss how to incorporate the 
displacement-distortion effect into our previously derived rate equations in terms of a joint 
electron-phonon representation using a generic quantum Langevin equation approach with 
a Markovian approximation. In this section, we also describe MacDonald's formulae for 
calculating the counting statistics, the first moment (i.e. current), and the zero-frequency 
second and third moments (i.e. the shot noise and the skewness, respectively). Then, we 
discuss in detail the counting statistics properties of a molecular QD for weak and strong 
distortion in Sec. III. Finally, a brief summary is provided in Sec. IV. 

II. MODEL AND THEORETICAL FORMULATION 
A. Model Hamiltonian 

We consider a simplified model of a molecular Hamiltonian constituted of two distinct 
parts, namely, electronic and nuclear parts within the Born-Oppenheimer approximation. 
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The electronic part is assumed to be constituted of one spinless level, Ed, with coupling to 
two electrodes [left (L) and right (R)] and also capacitively coupled to a gate electrode. 
Since an electron can tunnel into/out of this molecular QD from/into one of two electrodes 
in a transport process, there exist two possible electronic states: the molecular QD is either 
occupied by an electron (|1)) or is unoccupied (|0)). The nuclear part of the molecule 
is considered as a single harmonic oscillator. In this situation, there are two harmonic 
potentials corresponding to the two electronic states, which may differ from each other (we 
will term them the ground potential and the excited potential, respectively, in the following). 
The full Hamiltonian of the molecule can be written as 



H = Hieads + Hmol + Ht, (la) 



with 



Hieads = 2^£r;k4kC^k' W 



H^oi = ed\l){l\ + \^)(^\Hn, (Ic) 

N=0,1 

Hn = ^ + "^ujUx - d5mf, (Id) 
Ht = 5^(K,4|0)(1|+H.c.), (le) 

»7,k 

where cj^j^ (c,7k) is the creation (annihilation) operator of an electron with momentum k, 
and energy £:.^k in lead rj [r] = L,R). p (x) is the momentum (normal coordinate) of the 
molecular oscillator and m^, is the mass. When there is no additional electron occupying 
the QD, the frequency of this oscillator is ujq and its coordinate origin (the equilibrium 
position) is zero; but in the case of occupation, a distortion and a displacement occur, i.e., 
the frequency is changed to ui (which probably is different from ujq), and the equilibrium 
position has shifted through d, which is related to the coupling constant A = uj^d^frr^pj^J^ 
for the electron-oscillator interaction. is the tunneling amplitude between the QD and 
electrode 77. In this paper, we assume that is independent of the oscillator position, which 
is a good approximation for experimental realizations in most single molecule systems.-i^ii^i^ 
We use units with h = ks = t = rrio = 1 throughout the paper. 

Bearing in mind that the rate equation is written in the joint electron-phonon represen- 
tation, i.e., the direct product states of electron occupation number and phonon Fock state, 



it is convenient to rewrite the phononic part of the Hamiltonian, Eq. (fTal) . in quantized rep- 
resentation. We define phonon annihilation operators for the oscillator in two cases, = 
and 1, and phonon-bath, respectively, as: 

1 



/UN 



(2) 



In these terms, we have 

Hn = cuArajvOAT, (3a) 

The definition of Eq. clearly shows that there are two different Fock spaces spanned 
by the harmonic oscillator basis sets, \n)o and \n)i, corresponding to the ground and the 
excited harmonic potentials respectively, 

a^Naj^\n)N = n\n)N- (4) 

We will describe the electron-oscillator interaction in the tunneling Hamiltonian, Ht, using 
this representation. Consider a vacant initial electronic state of the molecule, |0), so its 
nuclear vibrational state is \n)o of the ground potential, i.e., the initial joint electron-phonon 
state (JEPS) is |0) ® |^)o (for notational convenience, we use 0„ to denote this JEPS). If 
an electronic tunneling event occurs, an additional electron is injected into the QD from 
one of the electrodes, so that the electronic state changes to In conjunction with this 
electronic transition, the nuclear vibrational state simultaneously changes to the Fock state 
of the excited potential due to external forces involved in the electronic transition in 
the molecule, leading to a final JEPS |1) ® \n')i (denoted by The total amplitude of 
the joint transition from the initial JEPS 0„ to the final JEPS !„/ (which is called vibronic 
transition in the literature) is equal to the product of both the amplitudes of electronic 
tunneling to lead rj, and the harmonic oscillator transition |?7,)o —>■ \n')i with the FC 
factor Fn'n, which is determined by the overlap integral of the wavefunctions of the involved 
harmonic oscillator states in coordinate representation, 

dxip*{x,uJo)ipri'ix - d,uji) = o(n|n')i, (5) 

where il)n{x,uj) denotes the eigenfunction of a harmonic oscillator with frequency uj in state 
In). Using Fan's IWOP method, we can derive an explicit analytical expression for this FC 
overlap integral as (Appendix)— 
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r 1 ( 2/5 ^inw+i)/2[nl^irm(l^y^' (-1 (n!n'!) 



1 + 13^ klk'l (n - 2k)\(n' - 2k')\ 

^ ^ k=0 k'=0 ^ ^ ^ ' 



2(72 + 7'2) 

/ I \ n-2k-n'+2k' , 

n' - 2k')\{-lT' ( 7\A^) i^:r^2'.7" ) , if n - 2A: > - 2fc'^ 



/ / N n'-2k'-n+2k / , n 



(6) 



2k<n' - 2k 



with P = a/co'i/cuo, 7 = i/tiJo'^ ~ (9 ~ ^/^o is the dimensionless electron- vibration 
couphng constant), and 7' = ^/uJld = jSj. L^{x) is the generahzed Laguerre polynomiaL 
Since the transport properties of the molecular QD depend crucially on the FC factors, it 
is essentially important in numerical calculations to compute these FC factors more accu- 
rately. The Eq. ([6]) provides a more convenient and precise way to accomplish this than 



the pure summation expression in Ref. 
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271 and the direct numerical integration in Ref. 
Besides, unlike the earlier paper,— where the FC factor is expressed as a product of two Her- 
mite polynomials, the Eq. ([6]) here is a more compact and tractable formula for numerical 



computations. Finally, the tunneling Hamiltonian Eq. (jlej) can be rewritten as 



Ht = J2J2(^vFnn'i^\0)\n)oi{n'\{l\ + Ii.c.). (7) 

rjjk n,n' 

B. Rate Equation 

Based on this representation of the Hamiltonian and following the same scheme we pre- 
viously proposed,^'' we have derived rate equations in terms of the JEPS representation 
of density matrix elements, Pqq = (|0) |n)oo(^| (0|) and = |n)ii(n| to describe 
unequilibrated vibration-assisted sequential tunneling. It should be noted that coherences 
between state with different phonon number can be safely neglected owing to the big energy 
difference between the two JEPSs 0(1)„ and 0(1)^ ifn^m and a small level broadening due 
to tunneling. In order to employ MacDonald's formula for calculating counting statistics of 
the tunneling current, we write the rate equations in a number-resolved version describing 
the number of completed tunneling events,— 

.n(0 _ \ ^ l-p- m{l) p_ 
Poo — 2-^ V L,nmP\\ R,nmHll 
m 

p+ n{l)l ( + _i_ ~ \ "(0 

nmPoO \ \^0,n ' ^0,n)PoO 
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. _ n+lfn . 4- n—l(l) /n \ 

+^o,n+iPoo +'^o,n-iPoo S (8a) 

.n(0 _ rp+ m{l) p+ m{l+l) 

PU — 2-^ V- L,mnHOO R,mnP00 
m 

_-p- "(0 _ /„+ , r-,- ^n"(') 
mnHll \^l,n "r '^l,nJPll 

+^l,n+lPll +^ln-iPn ^ (8b) 

in which the additional superscript (/) denotes the total number of electrons transmitted 
through the QD that arrive at the right lead during a fixed time interval t. Obviously, we 
have the relation p^j{t) = ^iP'jf\t) {j = 0,1). The electronic tunneling rates are defined 
as 

(9a) 

r- = Vr- 

nm / J ri,nm 

V 

= Y'^v'^nmi^- frii^d + muJi-nuJo)], (9b) 

V 

with = 27igri\Vrij\'^ denoting the tunneling strength between the QD and lead r] {grj is 
the density of states of lead rf). fr]{^) is the Fermi-distribution function of lead rj with 
temperature T, and the FC factors are 

The dissipation rates of the vibrational number states are 

^7V,n = 'CU.pnB{uJN){n + l), (10a) 

^Ar,n = G7p(?^i?(t^Af) + 1)^- (10b) 

Here, we assume the same vibration-bath coupling constant, tUp, for the vibrational modes. 
In the case of a pure displacement of the vibrational potential {(3 = 1), the FC factor becomes 
A.18I) . and thus the rate equations Eqs. (!8al) and (l8b1) reduce to our previous results in 



Eq. 
Ref. 



20. 



C. MacDonald's formula for couting statistics 

In principle, the full counting statistics of the tunneling current (all cumulants {{q^)) of the 
charge q{t) transmitted through a QD during a sampling time t) can be calculated employing 
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the number-resolved rate equations. However, the high moments are quite difficult to probe 
experimentally, due to the central limit theorem as well as their sensitivity to environmental 
influence.— On the other hand, as a classical current fluctuation, a Gaussian distribution 
has all its cumulants higher than the second one equal to zero, i.e. 

((g')) = for / > 2; while 

a Poisson distribution has all its cumulants equal to its mean (the flrst cumulant). The 
simplest measurement of the non-Gaussianity of the distribution q{t) is therefore the third 
current cumulant {{q'^)), which reflects the skewness of the distribution. Hence, we focus on 
the first moment {{q)) (the peak position of the distribution of transferred-electron-number), 
the second central moment ((g^)) (the peak- width of the distribution), and the third central 
moment in this paper. 

According to the definition, the three cumulants can be calculated as 

((g)) = 5^/P«(t), (11a) 
= 5]/^P«(^)-((g))^ (lib) 
= i2^'P^Ht)-m)){{q'))-{{q))\ (11c) 



where 



p^'\t) = j2[p-,!;\t)+pi'\t)] (lid) 



is the total probability of transferring / electrons into the right lead by time t, which satisfies 
the normalization relation YliP^^K^) — 1- 

In the long time limit, these cumulants ((g')) are proportional to the time, ((g')) = 
((X'))t + Q, with the current cumulants ((X')) and constant terms q (ci = 0). From the 
definition, we have^^i^i^i^ 

/ = ^A{q)) = (12a) 

CiZ t— >oo 

S2 = ^((g^))| ={{!')) 



dt 
d_ 
~dt 



I i— >oo 

(12b) 



Y^fp^'\t)-{tif 
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which are the average current and zero-frequency shot noise, respectively. Likewise, we can 
define the skewness as 

^3 = im\ ={{i')) 



dt 
d 
dt 



1 1— »oo 

(12c) 



l^P^^\t) - Ih'^ - 3^2 /t^ - 3/tC2 



t—too 



I 

To evaluate S2, Ci and 5*3, we define auxiliary functions G'^At) and H^At) (j = 0, 1) as 



H-it) = Y.^pfPii). (14) 

Their equations of motion are readily derived employing the number-resolved QREs, 
Eqs (EiD-dHbl), in matrix form: G{t) = MG{t)+gp{t) and H{t) = MH{t)+nG{t)+g' p{t) 
with G{t) = {Goo,Gnf, H{t) = {Hoo^H^Y, and p{t) = (poo,Pii)^ [here G,,- = 
{G%G]^r--Y, H,, = {H%,H]^r--f, and p,, = {p%p]^r--Y]. M, G, and H, Q' are 
easily obtained from Eqs. (I8al)-( l8bl) . Applying the Laplace transform to these equations 
yields 

G{s) = isI-M)-'gp{s), (15) 
H{s) = {sI-M)-^[nG{s) + g'p{s)], (16) 

with p(s) obtained by applying the Laplace transform to its constituent equations of motion 
using the initial condition p(0) = pst {pst denotes the stationary solution of the rate equa- 
tions). Due to the inherent long-time stability of the physical system under consideration, 
all real parts of nonzero poles of p(s), G{s), and H{s) are negative definite. Consequently, 
divergent terms arising from the partial fraction expansions of G{s) and H{s) as s — > 
entirely determine the large-t behavior of the auxiliary functions, i.e. the zero-frequency 
shot noise and the skewness, Eqs. ( ]12b[) and fll2cl) . 



III. RESULTS AND DISCUSSION 

In this section, we present the results of our numerical investigation of nuclear distortion 
effects on tunneling current, I [(in units of (F^ + Tii)/TLTfi], and on low-frequency counting 
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statistics, the normalized current noise (the Fano factor) F2 = S2/2I and the skewness 
-^3 = 5*3/2/ for the molecular QD. 

The vanishing of Wp = denotes no dissipation of the vibration to the environment, 
corresponding to maximal uneqilibrated phonon effects in resonant tunneling; whereas in- 
creasing dissipation strength, Wp > 0, describes the action of the dissipative environment 
as it begins to relax the excited vibration towards an equilibrium state (the later occurs in 
the hmit zup = 00). Throughout the paper we set cjo = 1 as the energy unit, /i^ = fiR = 
at equilibrium condition, and we assume that the bias voltage V is applied symmetrically, 
liL = —fJ'R = V/2 with symmetrical tunnel couplings, = F/j = F. 



A. Small shift of the equilibrium position 

To start, we consider cases involving a small displacement of the equilibrium position of 
the oscillator, A < 1. 

The case with an extremely small shift A <^ 1 (A = 0.01 in our calculations) corresponds 
to the situation in which charging the molecule induces only a distortion of the harmonic 
potential surface. It is already known that a pure displacement (no distortion, (3 = 1) retains 
the symmetry of the I-V curve with respect to the ground-state transition line {ea = 0), thus 
leading to a completely symmetrical ^d-dependent differential conductance, dl / dV .—^^^^ 
Moreover, there is no occurrence of negative differential conductance (NDC) for the single- 
level model with symmetric tunnel couplings.— 'i^'^ However, a charging-induced distortion 
{l3 7^ 1) will indeed break the symmetry of the excitation spectra for = (neutral state) 
and = 1 (anionic state), and thus also result in an asymmetric dl/dV-Sd curve due to 
the asymmetric FC factor -F„„/, as shown in Fig. [1], which plots dl/dV as a function of 
bias voltage V and gate voltage Ed for ui = 0.5 at T = 0.05. Surprisingly, this effect of 
pure distortion also gives rise to NDC even in the symmetric tunnel coupling model. These 
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results are consistent with the earlier prediction of Ref. 

In the following, we will provide detailed analyses of the effects of a pure distortion on the 
counting statistics of tunneling current and present an interpretation of the occurrence of 
NDC. We exhibit our calculated results without environmental dissipation, Wp = 0, for the 
I-V characteristic and the bias-dependent Fano factor, F2, for two systems, = —0.5 and 
0.5, with various different distortion ratios in Fig. [2] at temperature T = 0.02. It is evident 
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that (1) strong NDC occurs for the case of Ed = —0.5 and its pattern of NDC depends on the 
distortion ratio (3; (2) by way of comparison, weak NDC is found for the case of Ed = 0.5; (3) 
the distortion causes asymmetric behavior of the gate-voltage-dependent higher orders of 
counting statistics, the shot noise F2 and the skewness F3 (see Fig. 6 below); (4) distortion 
induces an enhancement of shot noise and even a giant Fano factor F2 up to 10^ for the case 
of Ed = —0.5. 

1. The case with < 

In an effort to understand these results, it is helpful to examine the bias voltage de- 
pendence of the joint electron-phonon occupation probabilities (JEPOPs) Pqq Pii- W^e 
plot the calculated results in Fig. [3] for the system having e^ = —0.5 with two different 
distortion ratios, ui = 0.7 [(a) and (b)] and ui = 0.6 [(c) and (d)]; and in Fig. H] for the 
system having Ed = 0.5 with ui = 0.7 as well. Consider first the results for Ed = —0.5 
and uJi = 0.7. It is clear from Fig. [3]^b) that the molecule is occupied by an electron (the 
state Iq) at the beginning and thus the current remains zero until V ~ I.Ocuq (since no 
additional electron can enter into the molecule from the left lead). When the bias voltage 
increases to l.Ocuo, the Fermi energy of the right lead is nearly equal to the energy of the 
JEPS lo, leading to a de-population of the state Iq and simultaneous populations of the 
states Oo and li. That is to say, in this situation, only three JEPSs Oq, Iq, and li are 
involved in tunneling events and all other JEPSs make no contribution to the tunneling, 
which is equivalent to a two-level QD with symmetric tunnel couplings {x = joo^ for the 
channel Iq and y = 7oir for the channel li) under a large bias voltage at zero temperature.— 
Therefore, we have Pqq ~ p^i — Pii — V^, the current is I = {x + y)T/3, and the Fano 
factor is F2 = [2{x^ + y^) + xy]/{9xy). Due to the FC selection rule for vanishing shift 
(x ~ 0.98r ^ y ~ 0), we obtain a huge Fano factor, F2 ~ 2x/9y ~ 10^. 

When the bias voltage increases further to IAujq, the JEPSs Oi and I3 become occu- 
pied (surprisingly), albeit their conventional resonant values should be = 2.00^0 and 
2(3u;i + Ed) = 3.2u;o, respectively. These unusual occupations can be understood quali- 
tatively in terms of vibration-induced cascaded single-electron transitions:— i^^i^S arbitrarily 
high vibrational excitations can in principle be accessed via cascades of single-electron tun- 
neling processes, if these processes become energetically allowed by applying an appropriate 
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bias voltage. For instance, if the Fermi energy of the left lead is located at O.Tcuo, the JEPS 
li is populated and, moreover, the single-electron transition, li —>■ Oi, is permitted because 
the bias voltage V = IAuq provides sufficient energy to activate it, Uq — (e^ + Ui) = O.Tujq. 
Concomitantly, the molecule is re-populated because the transition Oi I3 is also energeti- 
cally accessible {ed + Sui — cuq = O.Guq), and it is allowed by the FC selection rule as well. As 
a result, albeit the Fermi energy of the left lead is not aligned with the energy of the states 
Oi and I3, these states also become occupied [Fig. [3t^a,b)]. Another direct consequence of 
such cascade transitions is a decrease of the JEPOP p\i, and corresponding increase of the 
JEPOP pii- Except for the exclusive relation between the states Oq and Iq at the beginning, 
they attract each other at sufficiently large bias voltage due to the strong bidirectional cas- 
cade transition (Iq ^ Oq) rate 700 — 0.98. Therefore, the molecule is further de-populated 
at this value of bias voltage and thus the current increases [Fig. [21(a)]. Actually, it is evident 
from Figs. [3] and H] that p"^ = Pqq is always satisfied for all n at extremely large bias voltage 
since the permitted bidirectional transition 1„ ^ 0„ makes the two JEPSs act as a whole. 

Similarly, when bias voltage increases to F ~ l.Scuo (pl = O.QcJo), the JEPS I2 becomes 
occupied, leading to simultaneous decreases of p^i and Pqq, as indicated by the double- arrow 
in Fig. [3]^a,b). As a result, the molecule has a higher probability of being occupied by an 
electron, causing a suppression of current at V ~ I.Suq. 

The situation is a little different for the case with ui = 0.6 as shown in Fig. [3](c,d). For 
a small bias voltage V < I.Ouq, only three JEPSs Oq, Iq, and li contribute to tunneling 
as discussed above in the case of Ui = 0.7. We note the occurrence of the same values of 
current and huge Fano factor, F2 ~ 10^ [Fig.[2]^a,b)]. However, in contrast to the former case 
of ui = O.Tuq, when the bias voltage is increased to 1.4ci;o, there are still only these three 
JEPSs involved in tunneling, since the bias voltage 1.4^70 can not provide sufficient energy 
to enable the transition li — Oi [luq — {sd + uji) = 0.9ujo]. As a result, we find a further 
increase of p\i at this value of bias voltage, as indicated by the double-arrow in Fig. [3](c,d), 
and thus a corresponding decrease of Pqq leads to a NDC at the I-V curve. 

2. The case with Ed > 

The case of positive Ed is quite different. In Fig. HJ we plot the bias- voltage-dependent 
JEPOPs, Pqo and p^-^ for the system with Ed = 0.5 and ui = 0.7. It is evident that (1) 
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the QD begins to populate at ~ 2\ed\ — I.Ouq; (2) surprisingly, the states Oi and li are 
also nearly equally populated at this point (ppg — p\i because of strong cascade transitions 
between the two states), albeit that the Fermi energy of the left lead is not matched with 
the energies of the states Oi {I.Ouq) and li (1.2co'o); and there is not any cascade transition 
which can be activated to reach these two states. 

To explain these results analytically, we focus on the bias voltage region l.Ouo < V < 
2Aujq and we consider an approximative model with only four accessible JEPSs, Oo,i and 
lo,i, at zero temperature. The simplified rate equations without environmental dissipation 
read 

.0(0 _ p o(z) p o(z-i) p 1(0 p 1(^-1) 

Poo — ooPoo + ^ OOPll + J- OlPll + J- OlPll , 
Poo — ~U 10 + i iijPoo + lOPll + J- llPii 
Pii — -Uoo + J-iojPii +iooPoo + J- loPoo 

pI?) = -(2roi + rn)pS) + rnpS\ (17) 

with Vnm = '^Inm- In this case, with A = 0.01 and Ui = 0.7, we have Too = O.QSF and 
Tn = 0.95r > Toi ^ 1 X lO'^F, and T,^ ^ 7x lO-^F; and thus pgo ^ p?i ^ itlro//rll) - 0-4, 
Pqq ^ ~ 4j-]^_)_P^^^P^^-) ~ 0.1; the current is / ~ roo/2, and the Fano factor is F2 > 1. 
Moreover, we find a little monotonic decreases of poo in the large bias voltage region, which 
is responsible for the weak NDC of the I-V curve in Fig. [2]^c). 



3. Temperature and dissipation effects 

We have also examined the temperature and environmental dissipation dependences of 
the current and shot noise associated with a pure distortion effect. We plot the calculated 
results in Fig. [5] for the system with Ed = —0.5, A = 0.01, and uji = 0.7. It should be 
noted that the NDC and huge Fano factor are quite fragile in regard to the environmental 
dissipation of the vibrational mode: the NDC nearly disappears and the Fano factor becomes 
1/2 (the typical value for a single-level QD with symmetric tunnel couplings) with a weak 
dissipation rate Wp = 5 x IO^^cuq [Fig. [5](a,b)]; whereas the huge Fano factor is robust against 
increasing temperature (it remains 10^ up to a relatively high temperature over the range 
T = 0.2CU0 ~ O.ScJo). 
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4- Skweness 



We analyze the skewness or third moment of counting statistics of tunneling event through 
a molecule due to unequilibrated vibration and pure distortion effect. The skewness char- 
acterizes the degree of asymmetry of a distribution around its mean, i.e., the distribution 
of transfered electron number around the average current I here. A positive value of skew- 
ness signifies a electron transfer distribution with an asymmetric tail extending out towards 
more positive direction; while a negative value signifies a distribution whose tail extends out 
towards more negative direction. 

The skewness can be easily obtained, using MacDonald's formula as described above, for 
a single level QD with no coupling to a vibrational mode at a sufficiently large bias voltage, 



(r^-r^)^(ri + r|)+4rir| 
(r^ + r^)^ 



(18) 



37. It 



which is consistent with the result of the zero-frequency limitation reported in Ref. 
reads F| = 1/4 for an symmetric system, Tl = Tji, and Fg* = 1 if ^ {<^)Tji, which shows 
that the transport behaves as a uncorrelated Poissonian event in the strongly asymmetric 
case. For a molecular QD with coupling to a vibrational mode, we observe from Fig. [6] (we 
set A = 0.1 in these calculations here) that, due to distortion effect, the skewness exhibits 
(1) asymmetric properties with respect to e^; (2) a large negative value at the bias voltage 
region corresponding to appearance of NDC for the system with Ed < 0; (3) but a positive 
value for the case with > 0. For the case with = —0.5 and ui =0.7, the system reduces 
to an equivalent two-level QD with symmetric couplings, x and y (x ^ y), respectively, as 
indicated in section [Til A 11 at the bias voltage V < 1.4ti;o and zero temperature. Therefore, 
we can derive an analytical expression for the skewness as 

7x'^y'^(x + y) -2(x^ + y^) 7x 2x f x\'^ , , 

which gives a large negative value due to x ^ It is interesting to note that the emergence 
of a negative skewness in this system is completely stemming from the particular nonequili- 
brated population of the vibrational mode and the vibration-mediated tunneling rates at a 
certain bias voltage window, in comparison with that of a single level QD without coupling 
to a vibrational mode, in which the skewness is always positive at the whole bias voltage 
ranges. We believe that the appearance of large enhancement of the shot noise and the 
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negative skewness is a distinct signature of the effects of nonequilibrated phonon and vibra- 
tional distortion on the transport dynamics of a molecular QD. Likewise, we can also derive 
the skewness for the system with Ed = 0.5 and ui = 0.7 based on Eq. fll7p . the full form of 
which is too cumbersome to give here. Instead, we assume x = Tqq ^ Fn, y = Fqi ~ Fio 
and employ the fact, x ^ y, that leads to 

We also examine the role of environmental dissipation on skewness in Fig. [6](d). We 
observe that at equlibriated phonon situation ro^ — > oo, the skewness arrives at the typical 
value 1/4 of a symmetric system at large bias voltage region. Besides, the negative value of 
skewness for < is robust against increasing temperature (not shown here). 



B. Large shift of the equilibrium position 

1. Current and shot noise 

Here, we present numerical analyses of the effects on electronic tunneling of a large shift 
of the equilibrium position jointly with a finite distortion of the harmonic potential. It is 
already known that with no distortion, a strong shift of the potential can cause exponential 
suppression of the FC factors, and thus result in significant suppression of the vibration- 
modified electronic tunneling rates. This is to say that, a strong shift of the potential largely 
blocks electron tunneling in the low bias voltage region (FC blockade) and significantly 
enhances low-frequency shot noise.— The main effects of the distortion are the asymmetric 
dependence of differential conductance on Ed and appearance of weak NDC, as shown in the 
plot of the differential conductance as functions of Ed and bias voltage V for a system with 
A = 3.0 and ui = 0.8 in Figs. [7] and E^a). 

From Fig. [8]^a) and (b), it is evident that in absence of distortion {ui = 1.0), the large 
electron-phonon interaction strength induces significant current suppression and enhance- 
ment of shot noise due to the FC blockade.— >^»^'^ For instance, the effective tunneling rate 
between the states Oq and Iq is essentially reduced to Fqo = 1.23 x 10~^F owing to the FC 
factor for A = 3.0, which is responsible for electron trapping in the ground state of the QD, 
the state Iq, even at a resonant point of bias voltage, V > I.Ouq [as shown in Fig. M^a) 
and (b)]. It is also responsible for suppression of the current because II oc FooPoq 
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small bias voltage region. However, a distortion of the potential surface will increase this 
FC factor, Too- For example, the value of Too = 1-1 x lO^^F for uji = 0.6 and A = 3.0 is 
one order of magnitude higher than that with no distortion (see Table H]). This is the reason 
that we observe a relatively weaker electron trapping effect in the ground state [Fig. [n](c,d)] 
and an obviously enhanced current [Fig. [Hl(a)] for a system with a distorted potential in 
comparison with those having no distortion. 

Moreover, the distortion effect induces an oscillatory-type structure of the occupation 
probabilities of the JEFSs as shown in Fig. M^c) and (d), in comparison to the stepwise 
probabilities in the case of no distortion [Fig. [9](a) and (b)]. It is this oscillatory behavior 
of Poo Pii that causes the appearance of NDC in I-V curves as shown in Fig. [7] and 
Fig. [H^a) and (c). For the system with ea = —0.5, the molecule becomes empty at the 
bias voltage V = l.OcJo for both situations. Without distortion, the sequential transition 
lo ^ Oo ^ li ^ Oi ^ I2 ^ O2 ■ ■ ■ can be activated by the bias voltage V > I.Ouq 
and arbitrarily high vibrational excitations can, in principle, be accessed via the cascades. 
Nevertheless, only the probabilities of the four states Oo(i) and Oo(i) are exhibited in Fig. M^a) 
and (b) for the bias voltage region I.Oujq <V< S.Ocjq, while others are ignorable due to small 
transition rates. When the bias voltage increases to S.Oujq, some new cascaded transitions, 
e.g. lo ^ Oi ^ I3 ^ ■ ■ ■ , Oo I2 ^ O3 ^ ■ ■ ■ , are available and these additional 
channels induce stepwise increases of poo aiid the current. If an additional electron induces 
a distortion of the oscillator of the molecule (we set uJi = 0.6 here), the cascaded transition 
is quite different from that of the case of no distortion. When the bias voltage is about 
l.OuJo ~ 1.4co'o; the transition li ^ Oi is not allowed because the bias voltage is not strong 
enough to trigger this tunneling event, ujq — (e^ + ui) = 0.9ujq. Therefore, the infinite chain 
of the cascaded transition is broken, which means that in this bias voltage region, only 
three JEFSs Oq, Iq, and li are involved in transport. In this situation, the tunneling can 
be described by a simplified symmetric two-level model, as indicated in Sec. IIII A It with 
X = 7oo = 1-1 X 10"'^ and y = 701 = 9.5 x 10^'^ (see table [T]), which reveals a small current, 
/ ~ (x + y)T/3, an enhanced shot noise, F2 ~ 10, and a very small negative skewness [see 
Fig.in^a) below]. 

However, when the bias voltage increases above IAuq, another transition, Oq ^ I2, is 
activated with a stronger transition strength, 702 = 3.8 x 10~^, leading to a decrease of 
PQoand to population of the state I2, which also results in the opening of the transition. 
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l2 ^ Oi, with a stronger transition strength, 712 = 9.2 x 10~^ [because the energy change 
of this transition is quite small, uo — {sd + '2uJi) = O.Scuo], giving rise to population of the 
state Oi, as shown in Fig. M^c) and (d). This result is also different from that of the pure 
distortion case considered in Sec. IIII A ![ where the transition, I2 ^ Oi, is not allowed due 
to the FC selection rule and the state Oi is not populated until the bias voltage V = I.Suq is 
reached [Fig. [3]^c) and (d)]. If the bias voltage is increased to I.Suq, the transition, li Oi is 
opened, which causes (1) decrease of p\^; (2) further increase of Pqq; and (3) a feedback effect: 
increases of and pl^ due to the cascaded back-transition channel, Oi — I2 — > Oq — > lo- 
At the bias voltage V = 2.6ujq, two more transitions, Oq ^ I3 ^ O2 and Oq ^ I2 ^ O2, 
are activated, resulting in a decrease of Pqq again. At V = S.Ouq, the state Iq can directly 
transit to the state Oi: This new transition induces a strong decrease of p^i ^^id its strong 
feedback effect makes Pqq, pJq, and p\i all increase. To sum up, the newly opened cascaded 
transitions and their strong feedback effects are responsible for the oscillatory behavior of 
Poo and pn. 

Environmental dissipation will gradually destroy the NDC and the super-Poissonian shot 
noise, as shown in Fig. [H](c) and (d). Unlike the case of pure distortion [Fig. [S](a)], the 
current is obviously reduced by environmental dissipation. 

2. Skweness 

We have also examined the skewness of the vibrational-mediated electron tunneling in 
the case of a large shift of the equilibrium position of the oscillator, both with and without 
distortion. Figure [10] shows the skewness F3, in the case of no distortion, as a function of 
bias voltage for systems with = and A = 1.0 (a) and A = 3.0 (b). Unlike the case of pure 
distortion, there is no occurrence of negative values of F3. It is clear that in the equilibrated 
phonon situation (infinite environmental dissipation rate, Wp —>■ 00), the skewness has the 
typical value, = 1/4, for the symmetric single-level QD. For a moderate shift, A = 1.0, 
the skewness is slightly larger than in the small bias voltage region, but it is smaller than 
F| in the large bias voltage region. More interestingly, a finite dissipation rate obviously 
lowers the skewness below F| even in the small bias voltage region. Moreover, we find a 
significant enhancement (up to 10^) of the skewness in the case of a large shift, A = 3.0, 
which can also be ascribed to the FC blockade effect. 



17 



In the presence of distortion (Fig. [TTj) . the skewness exhibits asymmetry with respect 
to Ed and obviously gives rise to negative values in a moderate bias voltage region. For a 
system with Ed = —0.5 [Fig. [TTT a)]. distortion even results in a small negative skewness in a 
small bias voltage, as mentioned above. Environmental dissipation destroys these features 
of skewness (not shown here). 

IV. CONCLUSIONS 

In summary, we have analyzed, in great detail, the effect of distortion of the harmonic 
potential surface of a molecular oscillator (due to additional electron occupation) on unequi- 
librated vibration-assisted tunneling through a molecular QD in the sequential tunneling 
regime. This effect is modelled by a change of the vibrational frequency (phonon energy) in 
addition to a shift of the equilibrium position of the potential during electronic tunneling. 
To take this effect into account, we developed an explicit analytical expression for the vibra- 
tional overlap integral involving the harmonic wavefunction jointly with its shifted-distorted 
counterpart, which is known as the FC factor in literature. It is this factor that signifi- 
cantly modifies the tunneling rates during electron hopping between the electrodes and the 
molecule, and consequently it strongly influences the electronic tunneling properties of such 
molecular devices. In this paper, we employed Fan's IWOP technique to accomplish this 
task. 

Using these derived analytical expressions for the FC factors, we have established generic 
rate equations in terms of the JEPS representation to describe vibration-mediated tunneling, 
which facilitated examination of the roles of both shifting and distortion of the oscillator, as 
well as the roles of the unequilibrated phonon, and its dissipation to environment. Employing 
MacDonald's formula, we have calculated the current and its counting statistics, i.e., the 
low-frequency shot noise and skewness as well. 

Our analyses show that distortion has two main effects.— The first one is that, due 
to distortion-induced symmetry breaking between the excitation spectra of the two charge 
states of the molecule, i.e., the neutral state with = and the anionic state with = 1, 
the FC factors are asymmetric (7nm 7^ 7mn if n ^ m) and the molecular QD exhibits 
asymmetric tunneling properties with respect to > and e^ < 0. The second one is that, 
because the change of oscillator frequency can significantly change the cascaded transition 
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channels activated by an external bias voltage, it gives rise to different nonequilibrium 
vibrational distributions, and thus results in some new tunneling properties which are absent 
from the molecular QD without distortion. 

In particular, for the system of pure distortion (we suppose an extremely small shift 
of equilibrium position in numerical calculations), we have found that the molecular QD 
exhibits a strong NDC even at the condition of symmetric tunnel couplings, an enhanced 
shot noise with a huge Fano factor, and a large negative skewness in an anionic state; while 
it only has a weak NDC, a nearly Poissonian shot noise, and a large positive skewness 
for a neutral molecule. Two simplified models (and thus simplified rate equations) have 
been given to interpret qualitatively the reasons of different transport properties. In the 
presence of a large shift of the potential minima, the distortion also causes appearance of 
weak NDC and an enhanced shot noise, which can be ascribed to oscillatory behavior of 
the probabilities of the JEPSs due to the modified cascaded transition and strong feedback 
effects. The strongly enhanced skewness may be also negative at certain bias voltage. 
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Appendix: DERIVATION OF EQ. (ED 

For reference purposes, we briefly review the derivation of Eq. ([6]) following the Fan's 
IWOP method.-^ 

In a FC transition, e.g. electronic tunneling through a harmonic oscillator in the present 
paper, it is usually supposed that the initial state of the molecule (corresponding to the 
empty electronic state of the neutral molecule) is a standard harmonic oscillator. Ho 
[Eq. (]ld|) ]: whereas after an electronic transition occurs (the anionic molecule), the har- 
monic vibrational potential of the molecule. Hi, is distorted and displaced from that of the 
neutral one. Considering Eq ([2]), we have two different Hilbert spaces, |?7,)o and \n)i, and 
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the FC transition is described by Eq. ([5]). 

As suggested by Fan,— Eq. ([5]) can also be rewritten as 



^ nn' 



{x - d\n')Q 
o{n\ I dx -ojAx- d\n')o, 



(A.l) 

in which la;)^^^ and |x — (i)^^ are basis vectors in coordinate representation. Their expressions 
in Fock representation are 



^) exp ya;^ + V^xa^ - ^aj^^ |0) 



0; 



(A.2) 



with |0)o as the vacuum state of (ao|0)o = 0). It is convenient to define the dyadic 
(ket-bra) operator, S, as 



+ 00 



S = dx \x — d)u}^ 

J —oo 

X exp 



\ TT J 



oo 
t2' 



— {x - + V2^i{x - d)al - ^ 



yx^ + y/2uj^xao - ^ 



x|0)oo(0|exp 

Employing the operator formulae, 

|0)oo(0| = : exp(-ajao) :, 

exp(caQaQ) = : exp[(e'^ — l)alaQ\ : 



(A.4) 



(A.5) 
(A.6) 



: denotes normal ordering of the operators within the colons and c is a constant), we 



find 



S 



dx : exp 



X 



— {uJo + UJi) 



-yd^ - ^/2Ljldal - ^(ao + 4)^ 



2 12 

e 2(7^+7'^) gxp 



tanhx 1 2 sechx 



'-7% 



X exp 
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( c^n 1 sech x i \ 
X exp \ — tanh x + ^ 7 «o J , 



(A.7) 



with X = It is easy to verify that 1) S* is a unitary operator; 2) S\x)^^^ = \x — d)^^^, 

which means that S is an operator that transforms a state of the oscillator into a distorted- 
displaced state (FC transition); and 3) S'^aoS = qq cosh x ~ (4> ^i^i^ X ~ l' I = ol\- 
Using the operator S", the FC transition can be represented as follows, 



Fnn' = o(^|'S'^|'^')o, 



(A.8) 



which can be evaluated in the coherent state representation. Noting the relation between 
number states and the coherent state 

1 /dV 



Ho = ^|0)o 
V nl 



nl \da ^ ^ 



a=0 



we have 



Fnn' = (n!n'!)-i/2 



a \ I a \ 



da"* 



a*=T=0 



Considering two two coherent states a and r, we have 

ao\a) = a\a), (alaj = a*{a\, 
|a) = e"'^o|o)o, (a|r) = e"*^ 
{a\ : /(ao,4) ■ k) = 



(A.9) 



(A.IO) 



(A.ll) 



and, consequently, we can derive 



{a\S^T) = Aexp 



tanhx 2 2\ * 1 
— - — [a — T ) + a rsech x 



sechx, , * X 
H —['y a -it) 



where 



V2 



A = o(0|^t|0)o = sechi/2xexp 



2 '2 

7 7 



(A.12) 



(A.13) 



2(7^ + 7'^). 

Substituting Eq. ([02]) into Eq. flAnUD and using the definition of the double-variable 
Hermite polynomial,— 

Qm+n 



Qf^Qt'n 

min{m,n) 



exp[-U' + + <;t'] 



t=t'=Q 



m 



m—l^n—l 



1=0 



l\{m - l)\{n - l)\ 
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(A. 14) 



and the relation 



we find 



mini 



Ql+k 



rn—Ln—k 



F I 

J- Tj.r)/ 



2 /2 



2/3 



1 + 



(n+n'+l)/2 



[n/2] [n'/2] f l:±. 
A;=0 fe'=0 
xHn-2k,n'-2k' ^^'^^ 



(_l)fc'(^!^M)l/2 



(n - 2A;)!(n' - 2A;') 



/5 



.7 



(A.15) 



(A.16) 



Finally, employing the relation between Hermite polynomials and Laguerre polynomials, 



(A.17) 



we obtain Eq. ([6]). 

In particular, for /3 = 1 (no distortion), it is readily verified that Eq. ([6]) reduces to the 
well-known expression [see Ref. |20|, Eq. (9)] 



F I 

^ nn' 



(A.18) 



whereas for c? = (no shift of equilibrium position), this case is relevant to the theory of ac- 
tivity of non-totally symmetric vibrations in the electronic spectra of polyatomic molecules.— 
We easily obtain the Herzberg- Teller selection rule, F„„/ = if n — ra' 7^ even.— Moreover, 
we find 



nO 



(n!)i/^ 



exp 



2(72 + 7'') 

2 /2 \ "/2 



277' 



/ \ 1/2 



^2 _^ ^/2 



2(72 + 7'2) 



z/57' 



(A.19) 



and 



(n'!)i/2 



exp 



^2y2 



277^ 



1/2 
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in which Hn{x) is the single variable Hermite polynomial. These special results are identical 
to those of Englman's derivation.— 
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Table I: Several FC factors 'jnm (n, m = ~ 3) with A = 3.0 and uji = 0.6 (7nm / 7mn)- 





m = 


1 


2 


3 


n = 


1.1 X 10-3 


9.5 X 10-3 


3.8 X 10-2 


9.4 X 10-2 


1 


5.7 X 10-3 


3.5 X 10-2 


9.2 X 10-2 


0.132 


2 


1.6 X 10-2 


6.95 X 10-2 


0.114 


7.38 X 10-2 


3 


3.0 X 10-2 


9.6 X 10-2 


8.5 X 10-2 


1.0 X 10-2 
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Figure Caption 



FIG.l: Differential conductance dl/dV in the V-Sd plane for A = 0.01 and uji = 0.5 with 
Tl = Tn and with no environmental dissipation. The temperature is T = 0.05. The color 
bar gives the scale for the differential conductance. 

FIG. 2: (Color online). Distortion effect on the bias voltage dependent current, / (a,c) and 
Fano factor, F2 = 6*2/2/ (b,d) for the case of an extremely weak shift of the equilibrium 
position of the vibrational potential, A = 0.01 with Ed = —0.5 (a,b) and Ed = 0.5 (c,d), 
Tl = rii and with no environmental dissipation. The temperature is T = 0.02. 

FIG. 3: (Color online). Joint electron-phonon occupation probabilities, ppo (^''^) ^^"^ Pii 
(b,d), vs. bias voltage relevant to the parameters of Fig. [21 for (a,b) Ui = 0.7 and for (c,d) 
ui = 0.6. 

FIG. 4: (Color online). Joint Electron-phonon occupation probabilities, Pqq (a) and p"^ (b), 
vs. bias voltage relevant to the parameters in Fig. [2] for ui = 0.7 and Ed = 0.5. The thin 
lines in (a) and (b) denote the results for poo and pn, respectively. 

FIG. 5: (Color online). Calculated current I (a,c), Fano factor F2 (b,d), vs. bias voltage for 
the system with Ed = —0.5, A = 0.01, and uji = 0.7. (a,b) are for the results with various 
environmental dissipation rates Wp at the temperature T = 0.02; (c,d) are with different 
temperatures with no dissipation. 

FIG. 6: (Color online). Distortion effect on bias-voltage-dependent skewness, F3 = 5*3/2/, for 
the case of an extremely weak shift of the equilibrium position of the vibrational potential, 
A = 0.1 with Ed = —0.5 (a), (b), and 0.5 (c), and with no environmental dissipation, (d) 
Skewness vs bias voltage for Ed = —0.5 and ooi = 0.7 with various dissipation rates vup. The 
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temperature is T = 0.05. 



FIG. 7: (Color online). Differential conductance in the V-Sd plane for A — 3.0 and lui — 0.8 
with Fl = Fji and with no environmental dissipation. The temperature is T = 0.05. 

FIG. 8: (Color online). Current, I (a,c) and Fano factor, F2 = S2/2I (b,d) vs bias voltage 
for the case with A = 3.0 and a;i = 0.6. (a,b) are plotted for the cases with ea — —0.5 
(solid lines), 0.5 (dashed lines), and with no environmental dissipation. For comparison, we 
also plot the results without distortion (thin lines). (c,d) exhibit the roles of environmental 
dissipation for the system with — —0.5. The temperature is T = 0.02. 

FIG. 9: (Color online). Joint electron-phonon occupation probabilities, Pqq (^^^) ^^'^ Pii 
(b,d), vs. bias voltage; for (a,b) loi — 1.0 and for (c,d) loi — 0.6 with A = 3.0, Ed — —0.5 
and T = 0.02. 

FIG. 10: (Color online). Bias-voltage-dependent skewness, F3 = Ss/2I, for the case of no 
distortion of the vibrational potential with Sd — O and A = 1.0 (a), 3.0 (b), and with various 
environmental dissipation rates. The temperature is T = 0.02. 

FIG. 11: (Color online). Distortion effect on bias-voltage-dependent skewness, F3 = S3/2I, 
for the case of A = 3.0 with £d — —0.5 (a) and 0.5 (b), and with no environmental dissipation. 
The temperature is T = 0.05. 
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Figure 1: (Color online). Differential conductance dl/dV in the V-e^ plane for A = 0.01 and 
LVi = 0.5 with Tl = and with no environmental dissipation. The temperature is T = 0.05. The 
color bar gives the scale for the differential conductance. 
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Figure 2: (Color online). Distortion effect on the bias voltage dependent current, I (a,c) and Fano 
factor, F2 = 5*2/2/ (b,d) for the case of an extremely weak shift of the equilibrium position of the 
vibrational potential, A = 0.01 with = —0.5 (a,b) and = 0.5 (c,d), Ti = Tr and with no 
environmental dissipation. The temperature is T = 0.02. 
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Figure 3: (Color online). Joint electron-phonon occupation probabilities, Pqq (a,c) and (b,d), 
vs. bias voltage relevant to the parameters of Fig. [21 for (a,b) toi = 0.7 and for (c,d) toi = 0.6. 
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Figure 4: (Color online). Joint Electron-phonon occupation probabilities, Pqq (a) and p^i (b), 
bias voltage relevant to the parameters in Fig. [2] for uJi = 0.7 and £d = 0.5. The thin lines in 
and (b) denote the results for poo and pu, respectively. 
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Figure 5: (Color online). Calculated current / (a,c), Fano factor F2 (b,d), vs. bias voltage for the 
system with = —0.5, A = 0.01, and loi = 0.7. (a,b) are for the results with various environmental 
dissipation rates vjp at the temperature T = 0.02; (c,d) are with different temperatures with no 
dissipation. 
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Figure 6: (Color online). Distortion effect on bias-voltage-dependent skewness, = 3^/21, for 
the case of an extremely weak shift of the equilibrium position of the vibrational potential, A = 0.1 
with Ed = —0.5 (a), (b), and 0.5 (c), and with no environmental dissipation, (d) Skewness vs 
bias voltage for Sd = —0.5 and uji = 0.7 with various dissipation rates zup. The temperature is 
T = 0.05. 
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Figure 7: (Color online). Differential conductance in the V-ea plane for A = 3.0 and ui = 0.8 with 
Tl = Tji and with no environmental dissipation. The temperature is T = 0.05. 
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Figure 8: (Color online). Current, / (a,c) and Fano factor, F2 = S2/2I (b,d) vs bias voltage for 
the case with A = 3.0 and uji = 0.6. (a,b) are plotted for the cases with Ed = —0.5 (solid lines), 
0.5 (dashed lines), and with no environmental dissipation. For comparison, we also plot the results 
without distortion (thin lines). (c,d) exhibit the roles of environmental dissipation for the system 
with Ed = —0.5. The temperature is T = 0.02. 
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Figure 9: (Color online). Joint electron-phonon occupation probabilities, Pqq (a,c) and (b,d), 
vs. bias voltage; for (a,b) wi = 1.0 and for (c,d) uji = 0.6 with A = 3.0, = —0.5 and T = 0.02. 
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Figure 10: (Color online). Bias- voltage-dependent skewness, = S3/2I, for the case of no 
distortion of the vibrational potential with 8^ = and A = 1.0 (a), 3.0 (b), and with various 
environmental dissipation rates. The temperature is T = 0.02. 
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Figure 11: (Color online). Distortion effect on bias- voltage-dependent skewness, = 3^/21, for 
the case of A = 3.0 with = —0.5 (a) and 0.5 (b), and with no environmental dissipation. The 
temperature is T = 0.05. 
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